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Abstract 

We study the problem of dielectric response in the strong coupling regime of a charge transfer 
insulator. The frequency and wave number dependence of the dielectric function e(q,a>) and its 
inverse e _1 (q, u) is the main object of consideration. We show that the problem, in general, cannot 
be reduced to a calculation within the Hubbard model, which takes into account only a restricted 
number of electronic states near the Fermi energy. The contribution of the rest of the system to 
the longitudinal response (i.e. to e _1 (q, u>)) is essential for the whole frequency range. With the 
use of the spectral representation of the two-particle Green's function we show that the problem 
may be divided into two parts: into the contributions of the weakly correlated and the Hubbard 
subsystems. For the latter we propose an approach that starts from the correlated paramagnetic 
ground state with strong antiferromagnetic fluctuations. We obtain a set of coupled equations 
of motion for the two-particle Green's function that may be solved by means of the projection 
technique. The solution is expressed by a two particle basis that includes the excitonic states with 
electron and hole separated at various distances. We apply our method to the multiband Hubbard 
(Emery) model that describes layered cuprates. We show that strongly dispersive branches exist 
in the excitonic spectrum of the 'minimal' Emery model (1/Ud = U p = t pp = 0) and consider the 
dependence of the spectrum on finite oxygen hopping t pp and on-site repulsion U p . The relationship 
of our calculations to electron energy loss spectroscopy is discussed. 

PACS numbers: 71.35.-y, 71.27.+a, 79.20.Uv, 74.72.-h 
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I. INTRODUCTION 



The importance of many body effects for the description of dielectric response of in- 
sulating solids is generally accepted. In the neighborhood of the fundamental absorption 
threshold bound exciton lines and continuum excitons drastically change the spectrum. Re- 
cently it has been realized that the excitons in charge-transfer insulators (CTI) possess 
unusual features connected with the strongly correlated character of the antiferromagnetic 
(AFM) ground state in these compounds. The most pronounced peculiarity consists in the 
existence of exciton branches with substantially larger dispersion compared with the one- 
particle excitations.^ This behavior has a formal analogy with that of the Frenkel exciton^ 
that acquires a finite effective mass although both, electron and hole, have infinite masses, 4 
but in CTI it has a completely different origin. 

The experimental technique suitable for the observation of the exciton dispersion is the 
electron energy loss spectroscopy (EELS). 5 What is actually measured in transmission EELS 
experiments is the partial cross section^' 7 that may be decomposed into an amplitude factor 
and a dynamic structure factor 

The dynamic structure factor characterizes the linear response of the whole electronic system 
on longitudinal electric fields with the momentum q and frequency u ( the ionic contribution 
may be neglected for the considered frequency range of the order of several e V) . Pronounced 
peaks in <S'(q, a;) are called excitons.— They may correspond to discrete lines in the excitation 
spectrum of the solid or to resonances in the continuum part of the spectrum. 

For the theoretical description of excitonic features in conventional semiconductors and 
insulators the following scheme is used (The key references are Refs.yiQQQ). First, 
the quasiparticle excitation spectrum is found empirically 4 -**^ or from first principles . 1 ^ 11 
It is essential that the spectrum of the iV-electron system consists mainly of a continuum of 
electron-hole pairs whose electron or hole quasiparticle excitations are close to eigenstates 
of the (N + 1)- or (N — l)-electron system, respectively, with definite quasimomenta and 
energy— The quasiparticle spectrum is usually obtained from the self-consistent-field (SCF) 
approach. Next, the electron- hole interaction is taken into account. Then the problem for 
two quasi-particles interacting via the medium is solved. It is crucial that the ground state 
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may be viewed as an occupied valence band that is separated from the first excited state by 
an energy gap. 

For the CTI the above scheme should be revised beginning from the first step. The CTI 
has an odd number of electrons per formula unit. That is why the SCF calculations usually 
give a metallic ground state for the CTI and a gapless excitation spectrum. More elaborate 
SCF methods like the unrestricted Hartree-Fock method or its modern version LDA+U 
predict an AFM long-range order in the ground state and a gap in the excitation spectrum. 
Nevertheless, the nature of excitations and the ground state fluctuations are not caught 
by these approaches. Let us recall that any CTI remains insulating and shows excitonic 
features in optics and EELS spectra above the Neel temperature, i.e. in the paramagnetic 
state (with short-range AFM correlations). For strongly correlated systems the electron- 
electron interaction should be taken into account beyond the mean-field level. This is possible 
within the framework of Hubbard-like models in the restricted subspace of orbitals close to 
the Fermi level. 

The description of exciton physics in CTI is possible within the framework of one-band 
models , 14115116117 but more detailed and realistic information may only be obtained from the 
multiband Hubbard (Emery) model that explicitly includes ligand ion degrees of freedom. 

nnnn 

For quasi-one- dimensional cuprates the Emery model was considered in Refs. 18.19.20.21. 
The phenomenon of spin-charge separation that is characteristic to one dimension introduces 
a specific physics into the exciton formation. The mobility of a single electron or hole in 
ID is not suppressed by spin correlations and the exciton dispersion is comparable with the 
one-particle dispersion. 

In quasi-two-dimensional cuprates the situation is different. Due to AFM correlations the 
bandwidth of the one-particle motion is of the order of the AFM exchange integral J which 
is considerably smaller than the bare hopping t, and the exciton dispersion is of the order 



of t. The authors of Ref. 



221 proposed a qualitative physical explanation of the large exciton 



dispersion in layered cuprates: the propagation of an electron-hole pair does not disturb the 
AFM background, in contrast to the motion of a single electron or hole. Unfortunately, the 



calculations of Ref. 



22j do not support this idea since they give no dispersion in the absence 



of oxygen opjite Coulomb repulsion U r and direct 0-0 hopping t„, whieh was also pointed 
out in Ref. |18|. One should expect a qualitative description of the EELS spectr already 
for the 'minimal' version of the Emery model (1/Ud = U p = t pp = 0) which can be refined by 
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taking into account additional parameters like e.g. t pp or U p . On the contrary, the authors 



of Ref. |22j had to assume quite unrealistic parameter values to fit the experimental spectra. 

In this paper we outline an approach to calculate the longitudinal and the transverse di- 
electric response (Section II) for CTI within the framework of the multiband Hubbard model 
(Section III). Introducing an analogy to Wannier's excitonic representation^ we obtain a set 
of coupled equations of motion for the two-particle Green's function that may be solved by 
means of the projection technique^ (Section IV). It is substantial that the method allows a 
systematic improvement of approximations. In order to retain only the essential properties 
of CTI we first consider the 'minimal' version of the Emery model describing the Cu02 
plane of high-T c superconductors and their parent compounds. The model reflects the main 
features of CTI: the existence of two kinds of states - strongly correlated 'copper' states with 
the prohibition of double occupancy and uncorrelated 'oxygen' states; the undoped plane 
has one hole per unit cell, and is a quantum two-dimensional AFM insulator possessing no 
long-range order at any finite temperature; the charge excitation corresponds to the transfer 
of a hole from copper to the adjacent oxygen site. In Section V, we discuss the solution 
of our equations within the small exciton basis and its dependence on various parameters 
of the extended Emery model. We obtain an appreciable exciton dispersion and a rough, 
qualitative agreement with the experimental spectra already within the 'minimal' version of 
the Emery model. 

II. DENSITY-DENSITY CORRELATION FUNCTION 

The dynamic structure factor is related to the density- density correlation function 

S(q,co) ee — / dfe—* (n q (*)n_ q (0)> = ImA^), (1) 

2tti\ J-oo 7rexp — poj) — 1 



where 



= exp(-iqr) (<4 )S a r , s - (4 jS a r , s )) (2) 



q VN _ 

is the electronic density operator in the localized basis, the summation runs over all lattice 
sites r and orbital sorts s ; (...) means the thermodynamic average. For f3u > 1 we have 

S(q,u) « — — ImiV(q, u>), 

7T 
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where 

/*oo 

N(q,u) = «n q |n_ q » = -i / dfe"*" <[n q (t),n_ q (0)]> (3) 

Jo 

is the retarded Green's function that defines the inverse dielectric function 

Attp 2 

£ -l(q,u,) = l + — 7V(q,u;), (4) 
v c q 2 

with v c being the volume of the unit cell, and e is the electronic charge. We neglect here 
the local field effects (coupling with the Fourier components with q' = q + G, G being a 
reciprocal lattice vector). The function iV(q, uj) describes the response to the unscreened ex- 
ternal potential. It requests the account of the macroscopic electric field, i.e. the long-range 
part of the Coulomb interaction. The latter is responsible for the splitting into longitudinal 
and transverse excitons for small wave number which is analogous to the splitting into lon- 
gitudinal and transverse optical phonons.-^ The response to the total, screened potential is 
given by^ 

JV»(q,u;) = e{q,uj)N{q,u), (5) 

then 

Aire 2 

e(q,u;) = l jN s (q,u). (6) 

In the diagrammatic language the linear response to the total field may be expressed by 
the polarization operator where only irreducible graphs (which do not contain the contribu- 
tion of the macroscopic electric field) should be taken into account . 2 ^ 2 ? The random-phase 
approximation results in 



^(q,^) 

1 



ATI \ iv H4)"^ /«\ 
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which follows from (j3J) and (J3| and which is exact for q — > 0, as it was shown in Ref. 

The construction of Hubbard-like models for strongly correlated systems has an input 
from LDA band structure calculations where the screening of the long-range part of the 
Coulomb interaction is already taken into account. The Hubbard terms arise from the 
short-range residual interaction. Thus, the density response function N H (q, uj) calculated 
within the Hubbard model is an approximation to N s (q, tu).— In other words, it describes 
the motion of transverse (or 'mechanical' by Agranovich's 2 ^ terminology) excitons. 

Using the spectral representation we may write 

r°° T 1 1 Iuj' d.uj' f^o roo 

N s (q,z)= — ImiV s (q,u/) - -=/ +/ = N H (q, z) + N^q, z). (8) 

J o L vr \ z — uj Jo J ojo 



Here we bear in mind that the Hubbard model contributes to transitions in the low frequency 
region u < u>q with c^o of the order of the bandwidth, and the electrons of the rest of the 
solid are excited only at higher energies. In zero approximation we may assume that in 
the frequency region u > u> the electronic polarization of the rest of the solid follows 
immediately the external field 

JVootq, z) «iVoo(q,0). 

In other words, the Hubbard model is embedded into the medium with dielectric permeability 

, , 47re 2 



v c q 



-#oo(q,0). 



In fact, may has its own dispersion and may be quite anisotropic for layered or quasi- 
one-dimensional compound. In principle, it should be taken from e.g. LDA calculations (we 
have assumed that the rest of solid is uncorrelated) or from experiment. It is obvious that 
the peak positions of the loss function 



L(q, d) = — Im 



(9) 



and their intensity strongly depend on the value of £oo(q)- Usually one neglects the q- 
dependence and the anisotropy of £ ro , but it is a crude approximation, as well as another one 
which assumes e(q, 0) = const (see the discussion in Sec. V). For a quantitative description 
of EELS experiments the detailed knowledge of £00 (q) is necessary. Then the total dielectric 
function is 



e(q,o;) 



47re z 
v c q 2 



N H (q,uj) = e^Eu 



(10) 



and its inverse 

e -1 (q,o;) - c H 

For the dielectric function of the Hubbard model Eh the usual sum rule holds 



f- 1 f- 1 



u;Im£# (q, oj)dui 



uj\mE H {q,uj)duj = 

(, 2 EvoVcq 2 



h q ,H H 



71. 



where Hjj is the Hubbard model Hamiltonian and the operator n q acts in the subspace of 
orbitals which enter into Hh (i-e. the summation over s in Eq. (0) is restricted to these 
orbitals). Then for the total dielectric function approximated by Eq. (fTUj) we have 

ii 4ne 2 /it. a. 1 . i\ 1 



,71- 



E 2 JQ 



ulm.e(q,u))du>. (11) 



The factor arises due to the neglect of the frequency dependence of £ c 



III. MODEL HAMILTONIAN AND DENSITY OPERATOR 



As we have mentioned in the Introduction, we consider the 'minimal' Emery model that 
exhibits the essential properties of layered cuprates (l/Ud = U p = t pp = 0). Then the total 
Hamiltonian in hole notation reads 



H H = H + V, (12) 



where 



r >7 R.,o,7 

and where the Fermi operator p rrj annihilates a hole at site r of the oxygen sublattice with 
spin projection index 7, the Hubbard projection operator = Jr 7 (1 — n^) annihilates a 
hole with spin index 7 on a singly occupied copper site, where c?r 7 is the corresponding Fermi 
operator. The double occupancy of copper sites is thus excluded from (JT5J). H includes 
the on-site energies (A = e p — e<f, is taken as zero of energy), V is the p-d hybridization, 
a = x, —x, y, —y characterizes the direction of a nearest-neighbor vector a; the phase factors 
in V are absorbed into the definition of the operators p r (7 , , they do not change the 
exciton dispersion. 

Taking the limit L^/A — > 00 considerably simplifies the consideration and is a good 
approximation for weakly doped compounds in a wide range of values t/^/A > 2 (see e.g. 
Ref. Q. The conditions U p = t pp = are introduced for simplicity and may be easily 
relaxed (see Appendix), then the Hubbard on-site term for p-orbitals and the direct 0-0 
hopping of the form t pp = —t pp J2(i,j},-yPi itl Prj,-y is added to H . 

It is well known that the Hamiltonian ()12j) has an insulating ground state and is equivalent 
to the nearest-neighbor AFM Heisenberg model in the low-energy region. It means that 
charge fluctuations in V are strongly suppressed and that holes are localized. This fact 
becomes more apparent if we make a canonical transformation of operators of the form 

A eff = exp(-S)Aexp(S) = A+ A, S H , 

where 

R,o=±3:,±y,7 

Then Hh becomes 
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H, 



eff 



Ho - 4r z n + T E PR+a Ql , 7l PR+a It2 , 72 (ZrKv + Z 



R,7 R,al,a2,7 
•71 7O7] 

R+ga^R 



772,71 
J R 



R,o,7i 



(13) 



fsee also Ref. 



30 for the notation). Here p and Z mean transformed operators, J s is the 



AFM copper-copper superexchange interaction, g points to neighboring copper sites. Strictly 
speaking, the Hamiltonian (j!3)l is obtained under the condition t/A <C 1, and its parameters 
are r = t 2 / A and the AFM exchange J oc t 4 /A 3 . Nevertheless, it may be applied in a wider 
range t/A < 1 with renormalized values of r and J. 

The advantage of using the effective Hamiltonian (|13|) instead of the bare one ()12|) consists 
in excluding irrelevant zero-point charge fluctuations. Then the coupling of carriers with spin 
fluctuations that governs the low-energy physics of CTI becomes apparent. It is essential 
that the effective Hamiltonian ()13|) does not contain transitions between p- and Z-states, in 
other words, it never creates a particle-hole state out of the dielectric state. In this sense 
it resembles the starting Hamiltonian for the transverse exciton motion in conventional 
insulators. This allows to introduce an analog of Wannier's excitonic representation for the 
description of the electron-hole pair dynamics. 

The bare density operator 



£exp(-tqR) (2% - (2%)) + £«p(-iqs) 

R,cr S 



(14) 



(s runs over O sublattice) is transformed to 



N 



^exp(-zqR) 



R,<T 



2r+x E (pL^+^) 



a=±x,±y 
t 



+ 



(15) 



E eX P [-^( R + a )] PR+a^PR+a.a ~ ^ [(^R + ^R+2aJ PR+a Q ,a + k.C. 
H,a=+x,+y 

Note that in the second line of (|15|) (in the sum over the oxygen sublattice) a Q lies only in 
the same cell as R. Collecting the terms surrounding a Cu site, we have 

™ q = ™ q + -Lf^r E exp(-zqR) £ (Pn+* a ,a Z R + ^rW+b^) (1 - exp(-zqaj) 

= ^+ilE ex PN R ) E (^R,a a +^R,aJ (1 -exp(-zqa Q )) , (16) 
VJV A R o=dbx,±i/ 



where the operator 

7 

annihilates an electron-hole pair with minimal distance, and 
1 



(17) 



R. ( 



$>xp(-zqR) (Z™ - (Z£)) + ^exp(-zqs) [jfi - M 



;is) 



As we have mentioned above, the effective Hamiltonian (JT3J) conserves the number of particles 
in every band. Therefore n q = gives no contribution to Nh- Having the operator (|16|) we 
may proceed with the calculation of the density-density response function 



IV. ELECTRON-HOLE PAIR DYNAMICS 

The problem of the dielectric function (fTUJ) calculation is thus reduced to the calculation 
of the two-particle Green's function 

iV„(q, uj) = (($, + $t q |<t> q + $ \\ = (($ q |$t \\ + (($_ q |$t \\ (19) 



-ql^q 1 ^-l/Zu \\ t l | q. 

where 

$ q = tE ( X _ ex P(~ zc l a J) ^q,a, V'q.a = -7= EI exp(-2qR)?/> Ra , (20) 

and where we used again the conservation of particle numbers in the electron and the hole 
subsystem that excludes anomalous Green's functions like (($ q |$ q ))^. 
The equation of motion 

- (NO). = ( [*«> > + n*ff\ ( 21 ) 

generates more complex operators 

£n,a,p = EZ ^R+g Q -^R' ^ R +g«+ a ^'')'i > ( 22 ) 

71>72 



which annihilate states with an increasing separation between electron and hole, accompa- 
nied by spin fluctuations. The set of equations of motion will generate states corresponding 
to electrons and holes that are more and more separated and dressed by spin fluctuations. 
These states form a set similar to the excitonic representation for conventional insulators . 4 ^ 23 
The complication that arises in CTI consists in the strong interaction of both electron and 
hole with AFM fluctuations. The effect of this interaction leads to a strong renormalization 
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of the one-particle bandwidth, but it is partially canceled when electron and hole follow each 
other. 

The set of coupled equations (}2*T]) may be approximately solved by means of the projec- 
tion technique. We choose an operator basis B^i and the definition of the scalar product 
( Bq^B^j ). Within the operator subspace spanned by this basis we are looking for the 
approximate solution of the eigenvalue problem 

* q , H] = £ q * q , * q = £ Q(q)S qii . (23) 

i 

This leads, as usual for a non orthonormal basis, to the generalized eigenvalue problem 

5>(q)£u(q) = £5>(q)^(q) (24) 

i i 

where overlap and Liouvillean matrices 



S, 



'.i 



( 



Bi . B 



) = (B ti B}), L l3 = ( 



Bi,H 



B] 



) = ([B i ,H],B}) 



(25) 



depend only on spin-spin correlation functions for the system without electron-hole pairs, 
which is equivalent to the Heisenberg antiferromagnet. The correlation functions may thus 
be calculated from the Heisenberg model. 

The system (|24j) may be solved numerically. Then we find the eigenvectors 



*i = £<£(q)a«i,i, 



(26) 



where A is the number of the branch in the spectrum. In order to calculate iV#(q, u) within 
our basis, we should expand 



Then one obtains 



E k A (q; 



ql V q 



E k A (q; 



oj-e: 



with 



and finally 



g (q 



l x2 

A 



E cA (q)^« [1 - exp^qaj] 



47r 2 e 2 x . | A , 
v c q 2 Y 1 



u — E' 



-i 



(27) 



(28) 



(29) 



(30) 



Let us note that the projection technique allows to improve the chosen approximation step 
by step by enlarging the basis set. 
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V. RESULTS AND DISCUSSION 



We have restricted ourself to the minimal basis that describes the electron-hole pair with 
minimal distance. The basis contains operators (J 17)) and ()22|) with j3 = —a. Then the 
problem ()23|) has the dimension 8x8. The overlap and Liouvillean matrices are given in 
the Appendix. Spin-spin correlation functions were taken from the spherically symmetric 
treatment of an S=l/2 Heisenberg AFM model on the square lattice.— For a low temperature 
T = 0.1J and a vanishing frustration parameter p = 0.01 they have the following values 
(S R -S R+g ) = -0.33, (S R ■ S R+ga:+gy ) = 0.20, (S R -S R+2g ) = 0.17. For the on-site 
energy difference and p-d hopping we took the values A = 3.6 eV and t =1.3 eV, which are 
characteristic to all cuprates. 

Fig. 1 shows the dispersion of the imaginary part of the dielectric function £z(q,u) = 
Im£(q, uj). The oxygen on-site repulsion and the 0-0 hopping were neglected. We see 
strongly dispersive branches both in the [110] and in the [100] direction. 

As we have mentioned above, the comparison with the EELS experiment may have only 
qualitative character without a detailed knowledge of the background dielectric constant 
£oo(q). Fig. 2 shows the graphs for the loss function (multiplied by [^(O)] 2 in order to have 
approximately the same normalization as £2(q, w) according to (JHJ)) under the assumption 
that 22 

e(q,u = 0) ^e(q = 0,lu = 0) = 4.83. (31) 

Then the value of was obtained from Eq. ()30p. In this case the dispersion in the loss 
function reproduces essentially the dispersion in e(q,u). The peaks are slightly shifted to 
higher energies. In fact, the assumption (|3~T]) implies that the dispersion of Aoo(q, u) should 
follow the dispersion of A#(q, uj) in such a way that 

Ane 2 

e(q, u = 0) = 1 [N H {q, 0) + A 00 (q, 0)] = const 

as it follows from (|SJ) and (JEJ). 

In general, the interplay of A#(q, uj) and A 00 (q, a;) should be more complex. In order to 
demonstrate the strong dependence on the value of e^, we plot in Fig. 3 the same graph 
assuming a constant value £ ro = 2 for all q. We see a qualitative difference with Fig. 2 and 
we may conclude that the dependence of the loss function on is nontrivial. For quasi-one 
dimensional compounds a large values of ~ 8 was taken in Refs. Il^ol This means that 
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the rest of the solid strongly screens the long-range part of the Coulomb interaction between 
electrons that enter the Hubbard model. For this situation the poles of e(q, to) are very close 
to the poles of e -1 ^ u ) an d the shape of the loss function is close to the shape of £2(q? <^)- 
Note also that always e(q, 0) > as follows from Eq. (jl(Jj) . ()3Uj) and with = 8 one will 
receive unrealistically large e(q, 0). 

Let us now show some examples for the dependence of the dielectric function on various 
parameters of the model. For the reasons outlined above, from now on we consider only 
figures for £2(q? As it was mentioned in Section III, the model (|T2|) may be generalized in 
order to include a finite Hubbard repulsion U p on the oxygen site and a direct oxygen-oxygen 
hopping t pp . 

The dependence on t pp is not very strong for t pp < r. Let us recall that the addition of 
the 0-0 hopping to the Hamiltonian ()13j) within the parameter range t pp = (0.3 -j- 0.4)r is 
essential to describe correctly the angle-resolved photoemission spectroscopy (ARPES) of 
layered cuprates.— ^ Fig. 4 shows £2(q? ^) for t pp = 0.4r , and U p = 0. We see that the main 
difference to Fig. 1 consists in the redistribution of spectral weight between the different 
branches of the spectrum. 

The spectrum demonstrates a much stronger dependence on U p . Fig. 5 displays ^(q^) 
for U p = 4 eV. We have an almost vanishing dispersion of the lower branches. Let us pay 
attention to the fact that U p does not affect the single hole motion. Due to that reason its 
value is experimentally not well established. Our results show that the dielectric function 
dispersion is very sensitive to this parameter. 

Now let us discuss a very important peculiarity of our figures, namely the absolute position 
of the intensive peaks in the spectrum of ^(q, 1 ^)- At the T point (q = (0, 0)) we have one 
peak with the energy Er = A + 4r « 5.5 eV. Let us estimate the edge of the electron-hole 
continuum. It corresponds to the energy that is needed for the excitation of electron and 
hole which are independent on each other. The operator that annihilates such a state is 



Vi- h = e k+q /i_ 



i 12 eff 



(4 +q + ^ k )* q ""> 



where e k = — (— 4r — e k ~ J ) corresponds to the energy of a single 'electron' quasiparticle, that 
will be a complex spin-polaron corresponding to the coherent motion in the so called t — J 
model, which describes the Z subsystem in (113)) in the absence of holes (with t l ~ J = r). 
e k = A — e k is the energy of the coherent motion of the Zhang-Rice singlet, dressed by spin 
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fluctuations. The minimal energy that one needs to excite such a pair is 

E%£ ~ A - 4 in - (-At - ~ A - 2.2r « 2.6 eV . 

Here we have taken into account that the minimum of the spectrum in the t — J model is 
e min ~ ~~ 2r and that the Zhang-Rice singlet energy is e^ in ~ — 4.2r.— This result indicates 
that for layered cuprates the excitonic feature is immersed into the electron- hole continuum 
and represents a resonance rather than a discrete level. Of course, for the the final conclusion 
more detailed calculations should be performed within an enlarged basis set. 

Summarizing the calculated spectra we may state a rough, qualitative, agreement with the 
experimental curves^ already within the 'minimal' version of the Emery model and using 
a minimal basis set (Figs. 1-3). One may note a remarkable influence of the background 
dielectric function (Figs. 2 and 3). We found that additional parameters like e.g. t pp or 
Up act in different ways. That might improve a future parameter fit of the experimental 
curves. But the present accuracy is not sufficient for a reliable fit, which has to be reserved 
for the future. 

In this work we have not taken into account the intersite Coulomb repulsion which leads 
to electron hole attraction. This term has a twofold influence on the position of the excitonic 
feature. From one hand, it leads to an effective increase of the fundamental gap, from the 
other hand it contributes to the electron-hole binding. These tendencies are opposite to each 
other and should be thoroughly explored in a separate investigation. 
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APPENDIX 

Here we give the formulae for the matrices (}25|) for the Hamiltonian ( H"2|) that includes also 
a finite oxygen on-site repulsion U p and a direct 0-0 hopping t pp . The effective Hamiltonian 
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may be generally derived by the following procedure. Let be 



H = H + V, H = J2 \ m ) E m (ml, V = J2 \ n ) *nm ( 



ml 



(A.32) 



then the canonical transformation with the operator 

tnrn. 



s = 5» 



E m E n 



(m\ 



(A.33) 



gives 



Hq,S = — V , and finally up to the second order 



H e ff = exp(— S)H exp(S) = H + 



1 r 



v,s 



2 m \Em — Ej 



tnjtjm \ i i 

H — — {m\ . 

E n — Ej J 



(A.34) 



Within the operator basis containing (fTTjl and (f2"2"J) we have the equations of motion for 
the Hamiltonian containing finite U p , t pp 

[ipn,a,H] = E ^ n:a +Tj2^n,p+ru^n,a-a+[r u ^R +Sa -a,a-(r u -T)i> R+Sa _ a }-t pp J2 ^R,/3 

Here r M = t 2 /(A + C/ p ) and £o = A + 4r + 2(r — t u ) . It is convenient to introduce the 
notation 

UJ r = (^o' 7l ^g 1,72 '--^r 7r ' 7 )- 

7>7ii72v 

For r up to the third neighbors uj t is expressed via two point correlation functions 



— ~ + (Sr • S R+g ) , UJ Sa+Sfi - - + 2 (Sr • Sr. 



The Liouvillean and overlap matrices in k-space are then 



Sr ■ S R+ga+g/3 



(A.35) 



= 8 a pu g , (A.36) 

({[^ka, H],^}) = <W (r + r u cu 9 ) + (r - t pp ) (1 - <J Q/3 ) (1 - (A.37) 
+ ^,-/3[^ + exp(zkg a )(r u ^-r u + r)] , (A.38) 
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{{[^ka,-a,H],Jpl p }) = 5 a(3 {t u + TU g ) + (r - t pp ) u g (1 - 6 a p) (1 - 8 a -p) (A.39) 

+tf a ,-/3 {™<? + exp(tkg a )[r u - (t u - t)lu 9 ]} , 
({[£**,-<*, &],&(),-()}) = [<M r + T u u g ) + (1 - 6 afj ) (1 - S a _ p ) (r - tpp) uj g0+ga ] (A.40) 
+6 a -p {ru 2g + exp(ikg a )[T u u g - (r u - r)]} . 
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FIG. 1: The imaginary part of dielectric function £2(q, w) = Ime(q, ui) as a function of frequency 
ijj and wave vector q along two symmetry directions. For each direction the curve number n 
corresponds to q x = irn/lOa, a being the lattice constant, n = for the bottom curve. The 
parameters are A = 3.6 eV and t = 1.3 eV, oxygen on-site repulsion U p and 0-0 hopping t pp were 
neglected. 
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FIG. 2: The loss function L(q,w) = — Im [£ _1 (q, to)] (multiplied by [e^O)] 2 in order to have 
aproximately the same normalization as £2(q, ^)) for the same parameters as in Fig. 1 and under 
the assumption that e(q, lj = 0) « e(q = 0, uj = 0) = 4.83. 
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. 3: The loss function for constant = 2 for all q and the parameter set of Fig. 1. 
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FIG. 5: £2(q,o;) for t pp = and U p = 4eV. The other parameters are the same as in Fig. 1. 
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